Unconventional Josephson signatures of Majorana bound states 
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A junction between two topological superconductors containing a pair of Majorana fermions 
exhibits a 'fractional' Josephson effect, 47r periodic in the superconductors' phase difference. An 
additional fractional Josephson effect, however, arises when the Majoranas are spatially separated 
by a superconducting barrier. This new term gives rise to a set of Shapiro steps which are essentially 
absent without Majorana modes and therefore provides a unique signature for these exotic states. 



Majorana fermions comprise the simplest and likely 
most experimentally accessible non-Abelian anyon. An 
unambiguous demonstration of their non-Abclian ex- 
change statistics would be a great triumph for condensed 
matter physics, as this phenomenon reflects one of the 
most spectacular manifestations of emergence. Further- 
more, non-Abelian excitations provide the foundation be- 
hind topologically protected quantum computation [l|,y], 
with Majorana fermions playing a crucial role in pro- 
totype devices [3|-|7|. In the solid-state context, Ma- 
jorana modes were originally perceived as zero-energy 
states bound to vortices in p-wave superconductors 18[, 
and therefore are also associated with quasi-particlcs in 
the Moore-Read state M. More recent proposals employ 
topological insulators llCll-tl2f . half- metals in proximity 
to superconductors 131-1 15J , as well as spin-orbit-couplcd 
quantum wells [16|, ll7| and nanowires [18I421J to stabilize 
these elusive particles. Signatures of Majorana fermions 
appear in tunneling spectra and noise |22|,J23i, and more 
strikingly through interference effects [2J, |25[ . 

Josephson effects provide yet another important ex- 
perimental signature of Majorana fermions. Kitacv first 
predicted that a pair of Majoranas fused across a junction 
formed by two topological superconducting wires gener- 
ates a Josephson current [26 1 



I = tJm sin 
n 



(1) 



which exhibits a remarkable 47r periodicity in the su- 
perconducting phase difference (j)e ~ (f>r between the left 
and right wires. In stark contrast to ordinary Joseph- 
son currents, this contribution reflects tunneling of half 
of a Cooper pair across the junction. Such a 'fractional' 
Josephson effect was later established in other systems 



supporting Majorana modes [lO|, |ll|, [isl, |l9|, l27|, and 
in direct junctions between p-wave superconductors [28|. 
In this manuscript we demonstrate that two topological 
superconductors bridged by an ordinary superconductor 
with phase (j)m generically support a second kind of un- 



conventional Josephson effect with an associated current 



I = —Jz sin 
n 



(2) 



in the right or left superconductors, and twice that in the 
middle. This contribution arises solely from the fusion of 
spatially-separated Majoranas across the junction, and 
represents processes whereby a Cooper pair in the middle 
region splinters, with half entering the left and half enter- 
ing the right topological superconductor. We will derive 
this emergent term in Id Majorana-supporting systems, 
and propose several ways of measuring its effects. 

This novel Josephson coupling is derived most simply 
in a Id Kitacv chain. Consider a junction with Hamilto- 
nian H = Hi + Hr -\- SH, where the left/right supercon- 
ductors are described by p-wave-paired spinless fermions 
Ca,x {ct = £, r) hopping on an A^-site chain [26|, 



N-l 



Ha — — 



x=l 



{tcl ^Ca,x+1 + Ae^'^°'Ca,xCa,x+l + h.c). (3) 



j.^]^. The explicit 
in the Hamiltonians indicates 



Eq. (jSj adiabatically connects to realistic Majorana- 
supporting quantum wire Hamiltonians 18|, Il9l \2^ , and 
therefore describes their universal properties as well. Fol- 
lowing Kitaev, we express the spinless fermions in terms 
of two Majorana operators via Ca^x = h^~^~^ [Ib x + 
*7yi x)- When i = A, Eq. ^ maps onto a dimerized Ma- 
jorana chain: Ha = ~'itJ2x=i 1b,x1\ 
absence of 7^ j^ and 7^ ^ 
the presence of zero-energy Majorana modes localized at 
the ends of each superconductor in the junction. 

Let us now couple the two superconductors through 

&H ^ -im(c|jvCr,l + /l-C.)-A„(e"^'"C£,ArCr,lH-/l.C.), (4) 

where the two terms describe tunneling and Cooper pair- 
ing across the junction. These couplings combine the 
zero-energy Majorana modes residing at the junction into 
a finite-energy Andrecv bound state. Focusing on these 
zero-energy modes, one can write Cf^^v 



and Cr,\ -^ i\e 



'i4>^l1 



2^ IB,N 



7^ I , and define an ordinary fermion 



superconducting 
electrodes 




A, and (j) uniform the quasi-particlc gap is Eg 



\B 




FIG. 1: Topological insulator edge subjected to a magnetic 
field B and sandwiched by gates and superconducting elec- 
trodes. Majorana modes (red circles) localize at domain walls 
where the gap Egap = j \//i^ + A^ — B\ vanishes and the ar- 
gument in the absolute value changes sign. When the middle 
region is a trivial superconductor (S), and the sides form a 
topological phase (T), the novel Jz term in Eq. (O, with cur- 
rent ex sin(-^^^^ — 0m), accompanies the usual the fractional 
Josephson effect. This splits a Cooper pair in the middle elec- 
trode into two single electrons, injected via the two Majorana 
states in each topological segment. The same effect appears 
in spin-orbit-coupled wires in a T-S-T configuration. 



operator /''' = ^(7^ ^ + i"i\ {); SH then becomes 
SH ^ i2f\f-l){JMCOs[{cPe-<f>r)/2] 

+ Jz COs[((/)£ + (j)r)/2 - (f>,n]}- 



(5) 



to the fractional Josephson effect of Eq. ([T]) , while pairing 
A,„ across the junction produces the Josephson current 
in Eq. 1^ . Note that the sign of either current is dictated 
by the occupation number for the / fermion. and hence 
can be used as a readout method for qubit states encoded 
by the Majoranas [TTI.[29t. 

A more quantitative understanding is obtained by con- 
sidering more realistic models. Let us consider Majo- 
ranas localized on a topological insulator edge in proxim- 
ity to a superconductor and subjected to a magnetic field 
[lO|; a very similar analysis applies to quantum wires. 
In the Nambu spinor basis vp^ = {ip-f,i/ji,^],—'iljl), the 
Bogoliubov-de Gennes Hamiltonian for this system is 



n == vpa^T^ 



^r^ + A(cos0T^' -sin0r^) + So-^, (6) 



with V the edge-state velocity, p the momentum, B the 
Zeeman energy, and ct" and r° Pauli matrices acting in 
the spin and particle-hole sectors, respectively. We al- 
low the chemical potential p, pairing amplitude A, and 
superconducting phase (j), to vary spatially. 

Majorana states arise at interfaces between topologi- 
cal (T) and trivial (S) regions of the edge [10|. With fi, 



\/A^ + n^\ . When ^/A^ + ^^ > b the edge is gapped by 
proximity-induced superconductivity and forms a topo- 
logical phase closely related to that of Kitaev's model de- 
scribed above [lO|. In the trivial phase ■\/A^ + ji^ < B, 
and the magnetic field dominates the gap. We will study 
the T-S-T domain sequence of Fig. [l] which localizes 
Majoranas 71 at .t = and 72 at x = L. Each of the 
three regions, £, r, m, couples to a superconductor im- 
parting proximity strength A^/^/j, and phase (f>e/m/rj a-nd 
has a chemical potential fJ-i/m/r controlled by separate 
gates. (The main difference in the quantum wire case is 
that there creating the T-S-T domain structure needed 
to observe the unconventional Josephson effects discussed 
here requires the reversed criteria: ■\/A^ + p^ < B in the 
outer regions and a/A^ -|- p? > B in the middle region.) 
The Majorana-related Josephson effects result from 
hybridization between 71 and 72. When 71 and 72 are 
far apart {L — >■ 00), they constitute exact zero-energy 
inodes, and their wave functions decay exponentially in 
region s ^ £,m,r with two characteristic lengths: 



A 



s± 



\A,±^W^^, 



(7) 



with Ji\i =■ -f- and Jz = —f^ ■ Since the current in region ^ 
s is given by ^ i , ' , the fermion tunneling i,„ gives rise 



(we assume ps < B). For finite L, however, 71.2 
combine into a finite-energy state with creation oper- 
ator /^ = 5(71 + *72)- Roughly, each Majorana per- 
ceives the interface localizing the other Majorana as 
a perturbation, yielding a hybridization which is sup- 
pressed as a weighted sum of two decaying exponentials. 
This hybridization is again described by Eq. ([5]), with 



M/Z 



(j+e-^/^-+ ± J.e-i/^" 



An explicit cal- 



culation (see supp. material) for the symmetric setup, 
Pf ^ Pr = p, Ag = Ar = A > A,„ and pm = yields 



J4 



J- 



2A 



A(B+A)+M^ 



BA 



(8) 
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When Am = and the middle region is normal — which 
is the setup typically studied [lOl |l8[ — Jz = and hence 
only the Josephson term in Eq. ([l} appears. Turning on 
Am 7^ yields a nonzero Jz, and the second Joseph- 
son term in Eq. ([2]). Furthermore, since both Jz and 
Jm are dominated by the slowest decay length, they will 
generically be of the same order. For a quantitative 
estimate, consider the parameters pm — 0, pi^r = E, 
Am ~ E, Ai_r ~ V^E, B ~ 2E with energy scale 
E — O.lmeV. Assuming an edge velocity v = lO^m/s, 
for this choice we obtain Am-i- ~ 22nm, Am- ~ 66nm, 
and J± « 0.12meV. The effect then peaks at L « 50nm, 
which yields Jz « 0.022 meV and Iz = j-Jz ~ 5.3nA. 

These Josephson effects are simplest to understand 
conceptually when two additional Majoranas, 73,4, strad- 
dle the T segments of the edge as shown in Fig. [T] 
Let us define fermion operators /a = 2(71 ~^ *73) and 
fs ~ 5(72 + ^74), and assume that the corresponding 



occupation numbers are initially ua = 1 and ns = 0. 
Wc will further employ a 'pcrturbative' perspective and 
promote the superconducting phases to quantum op- 
erators conjugate to the Cooper pair number. One 
can then see that the Majorana operators in the term 
Ja,(2/V- l)exp(i^^) = zJM7i72exp(i^^) hop a 
single fermion across the S region, changing the state 
of the edge from (71^,715) = (1,0) to (0,1). At the 
same time, the exponential passes a charge e from side 
to side. The combination of these processes makes the 
term gauge invariant. The persistent superconducting 
current limit in this case is apparent when we consider 
an additional tunneling event which restores the parities 
of the T segments, moving a fermion back to the left but 
with a Cooper pair hopping to the right. A similar per- 
spective clarifies the role of the Jz term — the Majoranas 



in i J zl 1I2 exp 



.(. 



also change the parity 



of the two T segments, while the exponent removes a 
Cooper pair from the middle region and adds charge e to 
each T region (see Fig. [ij. 

Next, we discuss the crucial issue of measuring the 
new Josephson term in Eq. ©. The first and most di- 
rect possibility involves manipulating independently the 
phase differences 4>i — 4'm = *&l and (/)„i — (j)r = <&i?,, e.g., 
by inserting different fluxes in the two loops in Fig. [2ji 
(ignoring the voltage sources in the figure). By tuning 
$i = — $/j in a symmetric junction, one can probe the 
Jz Josephson term (driving current Jz sin $i on the mid- 
dle electrode) while canceling the Jm term. Such mea- 
surements, however, are highly challenging — they require 
careful flux control; the Majorana-relatcd Josephson cur- 
rent must be disentangled from the conventional 27r pe- 
riodic contributions; and the measurement must be con- 
cluded before the parity of the two Majoranas changes. 

A potentially more promising measurement scheme re- 
lies on Shapiro steps. In a regular Josephson junction, 
Shapiro steps arise from a combination of a dc voltage 
Vdc and an ac voltage Vac sin wi, which together generate 
a current I ~ Ij sin [(f>o + 2eVdct/h — {2eVac/fujj) cos wi]. 
Naively, this current averages to zero because of the con- 
stantly winding phase. This is not the case, however, 
when 2eVdc/fi — nco for some integer n — here a dc current 
component exists, producing a step in the V vs. / plot for 
the junction SOJ, |31|. For the fractional Josephson term 
in Eq. ([T]) , the Att periodicity leads to Shapiro steps when 
2eVdc/fi = 2nuj, corresponding to even Shapiro steps of 
a regular Josephson junction. The halved periodicity, if 
established, could provide a smoking-gun signature for 
Majorana modes. An inevitable conventional Josephson 
current, however, 'fills in' the missing ste ps, making it 
difficult to disentangle these contributions 28 1. 

The following three-leg Shapiro-step measurement cir- 
cumvents this problem and targets the Josephson term 
of Eq. ([2]). As shown in Fig. [2^, we envision a dc 
voltage applied to the left leg so that (fn — 2eVdct/fi, 
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FIG. 2: Three-leg Shapiro-step measurement scheme, (a) We 
envision applying a dc voltage Vdc to the left superconduct- 
ing electrode and an ac voltage with angular frequency uj in 
the left loop (which we model as an ac voltage applied to the 
middle electrode). A measurement of the dc current Ir in 
the right electrode will then reveal Shapiro steps stemming 
from the Majorana modes when the ac Josephson frequency 
2e Vdc/ft equals an even harmonic of lj. (b) Sketch of dIr/dVdc 
indicating the predicted Shapiro steps — note the crucial ab- 
sence of odd-harmonic peaks, which would appear in a con- 
ventional Shapiro-step measurement. 



while an ac voltage applied to the middle leg sets <j)m = 
— {2eVac/hjj)cosojt. Since the new Josephson term in- 
duces current in all three legs, a current measurement on 
the right lead will find Shapiro steps emerging only when 



2eVdc/fi ~ 2nuj 



(9) 



as illustrated in Fig.[5]D, without any odd-harmonic steps. 
This non-local measurement is insensitive to any parasitic 
two-phase Josephson terms, and therefore automatically 
eliminates most competing processes. Furthermore, it 
bears the advantage of being a fast dynamic measure- 
ment (since Josephson frequencies are typically in the 
GHz regime), which reduces its sensitivity to temporal 
fluctuations of the Majorana-state occupations. 

To verify the approximation methods used and to con- 
flrm the prominence of the Jz term in the three-leg 
Shapiro measurement, we also numerically analyzed the 
Josephson effects in a topological insulator edge. Fig- 
ure [3] shows that our analytical results [e.g., Eq. ([8])] in- 
deed agree very well with the exact numerical calculation. 
We also explored additional current contributions such 
as (5/sin((/)i + (pR — 20,„), which could obscure the Ma- 
jorana signature by producing unwanted odd-harmonic 
Shapiro steps. This term is independent of the Majorana 
modes, and can instead arise from conventional Bogoli- 
ubov states in the junction. In the limit of small pairing 
and tunneling over the middle segments, such a term re- 
flects a high-order process. Numerically, wc find that it 
is suppressed by at least an order of magnitude compared 
to the Majorana Jz contribution in the regime where Jz 
is substantial, i.e., when L is of order Xm±- 

By considering the full edge spectrum (including the 
Andrcev bound states and continuum states exactly) , we 
obtained the total Josephson energy of the domain con- 
figuration in Fig. [TJ 
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may also serve as a practical readout mechanism for qubit 
states encoded by these particles. 
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FIG. 3: Numerically determined coefficients of conven- 
tional Josephson couplings (Jl/r), Majorana-induced terms 
(Jm/z), and second harmonic of the Jz term (Jz,2). Our an- 
alytical estimates of J|J and J^ agree well with numerics. 
The energy unit is E and the length unit is ^ = v/E. The 
parameters are fii^r = E,fj,m = 0, Ai,r = VSE, Am ~ E, 
and Bi^r = Bm = 2_E. The characteristic lengths are 
\m+ = C/3 and Am- = S,- For E — O.lmeV and v = 10*m/s, 
the length unit is ^ = 66nm and the maximum current is 
Iz = j^Jz ~ 5.3nA. 



+ Jm C0s[((/)^ - 0r)/2] + Jz COs[{(j)t + (j)r)/2 - (f)^] 

oo 

+ ^ Jz,«cos[n((0^-h0,,)/2-0,„)] + --- (10) 

Here Jl/r are conventional Josephson terms (to which 
the three- leg measurement is insensitive), J^z/z are the 
Majorana-induced contributions, and Jz,n denote the 
(unwanted) higher harmonics of the Jz term. As Fig. 
[3] illustrates, Jm dominates for L <C Am+. while for 
Am- ^ L > Xm+ the Jz term becomes comparable, 
enabling the three- leg Shapiro-step measurement. The 
higher harmonics Jz,n are at least an order of magnitude 
smaller than Jz in this regime and can be neglected. 
For L 3> Am- the Majorana signatures are strongly sup- 
pressed as expected. 

In this manuscript, we explored a new Josephson ef- 
fect that arises when a pair of Majorana fermions fuse 
across a junction formed by two topological supercon- 
ductors separated by an ordinary superconductor. The 
Majoranas in this setup enable Cooper pairs injected into 
the barrier superconductor to 'splinter' into the left and 
right legs of the junction — a process which would ordi- 
narily be prohibited at low energies. While Majorana 
modes can also give rise to a novel fractional Josephson 
effect in T-normal-T junctions, we argued that an im- 
portant advantage of our setup is that here one can more 
readily isolate the Majorana-mediated Josephson current 
through Shapiro-step measurements. The experiments 
we proposed could provide a relatively simple and un- 
ambiguous detection scheme for Majorana fermions, and 
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SUPPLEMENTARY MATERIAL 
PERTURBATIVE CALCULATION 

Let us pursue here a detailed calculation of the Joseph- 
son coupling across the Majorana junction described in 
Fig. [TJ We will first find the wave functions of the Ma- 
jorana states localized on each domain wall, ignoring the 
existence of the other interface. We will denote these 
states as \L) and \R). Next, we follow the usual proce- 
dure for finding tight-binding states and Hamiltonians. 
We first calculate the overlap matrix, M^p ~ (a| /?) with 
a, /3 = L, R, and the Hamiltonian matrix within this 
subspace, ha/3 = (a| 7^1/3). It is easy to see that the 
approximate hybridization Hamiltonian is then given by 



H„^aJ =- M-^'^hM-^'\ 



(11) 



A(x) . 


= A^e(a;^ 


-L) + A^(l-e(a; 


-L)) 


^l{x) -- 


= ^rB(x- 


' L) + ij,^{l - Q{x ^ 


-L)) 



Single Majorana solution at a; = L. Wc first solve 
for the zero-energy cigcnstatcs of the Hamiltonian ^ 
with parameters: 



(12) 



The solution has the same form on the two sides of the 
domain wall, but with different parameters. We denote 
the side of the domain with the index s being s = r, m for 
right or middle. By squaring the Hamiltonian and look- 
ing for momentum values yielding zero energy states, we 
find two imaginary momenta on each side, which corre- 
spond to the spatial decay constants AJj. given by Eq. 
^ . The wave function associated with each side of the 
domain is: 



\r) = ^six) = Rs+i^s+e ^»+ + Rs-i's-e 



(13) 



with Rs± being four complex numbers determining the 
amplitude of the wave functions corresponding to the two 
decay lengths, and with ^s± being four, four-dimensional 
vectors, which when 05 = are given by: 



/ 



^, 



rnib 



1 



\ 



±i 



\T^ 



-»C,r 



/ 



i'rd 



1 

i 



\ 



(14) 



s+»^/^^ 



with cxp(iCs) = — — g ^, for s = £, r, m. These 

solutions are the building blocks for each Majorana state. 
In order to obtain what the wave function becomes when 
(ps (the phases of the superconducting electrodes) deviate 

from zero, we can apply the rotations: [/^ = exp ( i^- 
such that: 






UcbA 



S±- 



(15) 



Obtaining the Majorana solution follows from matching 
the boundary condition of the solutions, and from them 
finding the coefficients Rs±. 

To avoid the complicated expression that could arise in 
the most general case of Majorana coupling, we concen- 
trate on the case where A^ = A^ > A„i and ^Lj. = ng, = jj, 
and ^m = 0. This choice does not constitute a substan- 
tial loss of generality, and is useful for grasping the results 
of our calculations. A straightforward but rather tedious 
calculation leads to the following solution for the ampli- 
tudes of the decaying waves of the right Majorana state 
under the above assumptions: 



_ 2sinC,- 



Rr+ 
Rr — 



I sm 



cos 

1 



(^-^) 
(^ ■ 



(16) 



By symmetry, we can infer the structure of the left 
Majorana, which is localized about x = 0: 



\L) = -^fXx) = Ls+i's+e ^'■+ + Ls-^s-e' 



(17) 



The amplitudes Lg± also depend on the phases on the 
left and middle segment of the wire in a similar way: 



L 



■?n+ 



2 sin (^£ 



-I sm 



0<!-0TT 



COS 

l-ie~''^g 

1 



(^ 



(18) 



From the above results, and under the symmetric 
choice of parameters, we can compute the overlap ma- 
trix. Map = (a| /3). Neglecting exponentially suppressed 
corrections, we obtain the following form: 



M, 



ap 



vS, 



-B+A„ cos(0o 



"/^ [ 2(B2-A^) 



^) 



2A^(A2,--B2+^2 



(19) 



with V being the spin-orbit velocity. 

The coupling between the Majoranas could be calcu- 
lated perturbativcly by considering the two domain walls 
juxtaposed. For instance, while the left Majorana is an 
exact zero-energy eigenstate of the Hamiltonian 



the existence of the right segment of the wire perturbs 
this wave function, with the perturbation potential being 

Vr^ iHr-Hrn)0{x-L). 

Similarly, we can write H = T-Ln + Ve with Vi = {Hi — 
T~l,m)9{~x). This perturbation induces a hybridization 
matrix between the left Majorana and the right Majo- 
rana: 



h^ 



In our case, 



{L\Vr\L)^Q {L\Vr\R) 

{R\Vr\L) {R\Vi\R)=Q 



(20) 



the spin and particle-hole sectors, respectively. We al- 
low the chemical potential /i, pairing amplitude A, and 
superconducting phase (/), to vary spatially. In region s 
(with s = l,m, r), the parameters (/x. A, </>) = (/i^, Ag, (j)s) 
are constant. Without loss of generality, we assume 
(/),„ = to be a reference of superconducting phase. 

The Josephson effects in the TST junction has both 
bound states and continuum contributions. In the fol- 
lowing, we first present the procedure to compute the 
exact interaction energy E between two Majoranas, and 
then provide the formalism to calculate the energy con- 
tribution from the continuum. 



Vr — [{Ar COS (f)r — A,n COS (f>m) T^ — 

(Ar sin0r — Am sin 0m) r'^ — /^r'^'z] ^{x — L) 
The perturbation matrix we obtain is: 



-i/A™ 



e--^/^+ sin ■ 



COS '^'' n cos ■ 



e-0n 



(21) 



with V an unimportant phase. 

We arrive at the final answer for the Josephson cou- 
pling using Eq. pip . The result indeed coincides with 
Eq. (0): 



n 



JM 



(2/V-l) 



J+e-^/^+ sin ^^^^^ sin '■ 



-J_e-^/>^-co&i^^^^cos- 



= (2/t/-l)(jMCOs^^ 
with the constants J± being: 



Jz cos 



2 



where Mr 



B 



Mr 



(22) 



(23) 



is the av- 



2A,-(A2-B2+^2 

erage of the overlap matrix [Eq. ((T9)) ] diagonal elements, 
dropping the cosine term. The cosine term in the over- 
lap will produce additional harmonics of the Majorana- 
Joscphson term but will not qualitatively change the an- 
swer we obtained. The J± terms give rise to to the pre- 
viously explored Majorana- Josephson term, Eq. ([T]) and 
to the new zipper term, Eq. ([2]). 



NUMERICAL CALCULATION 

We now detail the procedure of our numerical calcula- 
tion. In the Nambu spinor basis ^^ = (-0-1-, -0_[_, ■)/;!, —^'I), 
the Bogoliubov-de Gennes Hamiltonian for this system is 

n = vpa'-T'- -^iT^'+A (cos (t>T^ - sin (/)t^) + Bct^, (24) 

with V the edge-state velocity, p the momentum, B the 
Zeeman energy, and cr" and r" Pauli matrices acting in 



Bound state energy 

For TST configuration, there are two Majoranas at 
interfaces between topological and trivial regions. The 
finite separation leads to a finite interaction energy 
E = Eint between these two Majoranas, with spatial- 
dependent wave function ^ = ^ (cc) satisfying the equa- 
tion 



•H* = E^^. 



(25) 



We will solve the interaction energy E — Eint by match- 
ing the boundary condition of the wave function. 

First, we replace the momentum operator p with —i-§^, 
and obtain the linear differential equation associated with 
energy E 



^m{x) = GE^{x), 



(26) 



with 4x4 matrix 



E 



Ge = i—<T~A — a'- (coS(/)T^ + sin^T^^) 

V V V V 

(27) 
In region s (with s = l,m,r), the parameters (^, A, (p) = 
{fig, As, (ps) E^re constant, and the matrix G^ has eigen- 
system 

G^^^^-^' = nfu^^ (28) 

with eigenvalues k,^ and eigenvectors u,^ for j = 
1, • • • ,4 and s = l,m,r. 

Then, we expand the four-component wave function 
^ (x) in terms of eigenvectors u-^ . We arc interested in 



the localized state with E < E, 



(':'■) 



In the left region, 



there are two localized modes (ReK-j^ 2 > 0) ^^^ the two 
divergent modes (RcKg 4 < 0). Similarly, in the right re- 
gion, there are two localized modes (ReKj^ 2 < 0) and the 

two divergent modes (RcKg 4 > 0). The wave function 
with two localized Majoranas consists of localized modes 

J2 cfe^^P^'ifl^ for2;<0 

'^''>)e4'-'(-^)4'-)forx>L ■ ^^^^ 

J = l,2 



*(a;) 




FIG. 4: The function det Mb vanishes when E = ±Eint, 
which can be used to numerically find the interaction energy 
between the Majoranas. 



In order to match the coefficients associated with left 
and right regions, we integrate the wavefunction over the 
middle region and obtain the condition 



Y, 4''4^ = ^{L)^ 6°-"'% (0) = Yl 4"^''" "^ 






which can be written as 



J = l,2 



M, 



„('■) 
V c^ J 








(30) 



(31) 



with 4x4 matrix 

Me- 



e^^-'^^'M (eG^"'i4') 



4') (-to ■ 

(32) 
The necessary condition for non-zero solution is 



det Me = 0, 



which can be used to numerically determine the interac- 
tion energy Eint ■ As illustrated in Fig. |4l the function 
det Mb vanishes at E ^ ztEint- (There is a technical 

(I r) 

subtlety associated with the fact that G^' is not a Her- 
mitian matrix. For some fixed values of E, the eigenval- 

(l r) 

ues of G^' have multiplicity larger than one, and the 
eigenvector u- ' might be a zero vector, which may also 
lead to spurious solutions with vanishing detMg. This 
issue can be resolved by using a polynomial discriminant 
to identify and remove these spurious solutions.) 



Continuum contribution 



the scattering matrix S^ ~ SE{4>i,4'r), which can be 
computed by matching the boundary conditions for all 
incoming and outgoing modes. Once we know the scat- 
tering matrix, we can use the Fumi's sum rule to compute 



the continuum contribution to the system energy [32l . |^ 



rlE 

WU,.(l}r)= I — In [det [S£(0/, </),)]]. (34) 

■p 

The continuum contribution consists of many Fourier 
components 



Wi^l,(Pr) 



with W„ 
terms are J 



E 

ni .nj- — — cxD 



Wni,nrCOs{ni4:i+nr4>r) (35) 



W-m.-rw Then conventional Josephson 
L/R = W^i,o and Wo,i, and the even harmonics 
of the zipper terms are Jz,2n = Wn,n- In the following, 
we provide the formalism to compute the scattering ma- 
trix Se i(f>h(f>r) ■ 



A',r) 



For energy E > Egap , there are the propagating modes 



(RCK 



O.r) 



0), with momentum p 



{l,r) 



= ImAcJ- ''^ . Sup- 
pose there are four incoming modes ( wj^ , Uj , u]^ , ^2 
with p\2 > ^ and p-[ 2 < 0> and four outgoing modes 



-.(/) ^(/) ^(r) ^(r) 

ig ,1(4 , Ug , U4 



with pf^ < and p^'] > 0. The 
wave function can be written as a linear combination of 
all these modes 



E 



^e^'-v!j^ 



for X < Q 



Y, cy -'e^ '^ ^'ity for X > L 



j=i, 



In order to match the coefficients associated with left and 
right regions, we integrate the wavefunction over the mid- 



(33) die region and obtain the condition $ (L) = e 



G'. 



^*(0). 



The relation between the amplitudes of incoming and 
outgoing modes is 



/ 4'^ \ 



.(') 



ME,,n 

with 4x4 matrices 



(r) 

1 

,(r) 



= M 



E.out 



( 4'^ \ 

(r) 



(37) 



M 



E,in 



.G'. 



^lAI) 



''lJI) 



-1i< 



-rir) 



(38) 



We now consider the energy contribution from the con- 
tinuum. The continuum states can be characterized by 



M 



E.out 



=.G' 



'i„-rW 



''Lji) 



Jr) 



r('-; 



(39) 



The scattering relation is 
/ 






Pi 






\ 



1/2 



(r) 



(pr 

with scattering matrix 



-P3 

-pi' 



(0^^^^ JO \ 
'-a 



'-4 



-prV''cr) 



VH^')'''4^^ 



(40) 



ss=pi/^M-^„„,.M^,„.p-^/^ 



(41) 



where P,„ = Diae 



(0 (0 (»-) (»■) 

Pi ,P2 ,Pl >-P2 



and P„ 



(0 (l) {r) (r) 
Ps ,P\ ,P3 ,P\ 



-Diag 

vation of current is 



The requirement of conser- 



E^i 



(0 



,(0 



Ep. 



Sr) 



(42) 



which ensures the unitarity of the scattering matrix 



^E^E 



(43) 



Hence, dot [S^] = c^^-Jb ^^^ ^^ hi [det [Sij (0;, c/)^)]] = 
—Se {4>ii4'r)- NumericaUy, we just need to compute the 



quantity Se {4'1,4'r) and the integral 



^°° dE 
W (</.,, 0,)= / —5E{4^u4^r)- (44) 



The continuum contribution W {(j)i^(j)r) has 27r period- 
icity in both (pi and 4>r, with Fourier decomposition of 

W {(j)l,(j)r) = \ I]^,n, = -oo W^"i,n, COs{ni(j)t + nr4>r), with 

Fourier coefHcients of Wm^n,.- The relevant Fourier com- 
ponents are Jl/r = Wi^ = Wo,i, and Jz.2n = W^.n for 
n = l,2,---. 

There is one subtle issue in the computation of the 
scattering matrix. There are four propagating modes 



for E > 



are two propagating modes and two localized modes 
for B('''-' + A/(M('.'-))' + fA('.-))' > E > E. 



but there 



il,r) 



gap 



B(M-) 



(/xC^-) 



(AC-'-) 



In the latter case, we 



need to compute the effective scattering matrix that are 
projected to the subspace spanned by the propagating 
modes. 



